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PREFACE 

This  report  deals  with  one  aspect  of  RAID's  continuing  investigation 
of  the  effects  of  the  space  environment  on  artificial  earth  satellites. 

The  problem  dealt  with  here  is  the  naturally  occurring  charge-up  or 
electric  potential  attained  by  a  satellite  as  it  passes  through  the  ionized 
regions  above  the  earth.  The  potential  may  have  important  effects  on  the 
response  of  instruments  in  space  probes,  on  radiations  inside  manned  space 
vehicles,  and  on  satellite  lifetimes.  Thus,  the  present  report,  while  of 
only  theoretical  interest  in  Itself,  should  ultimately  have  application 
in  these  areeis  of  space  studies. 


SUMMARY 


As  an  aid  in  the  detennination  of  the  electric  potential  of  natvirally 
charged  satellites,  the  concept  of  the  idealized  sheath  introduced  hy 
Langmuir  and  Mott-Smith  is  studied  in  some  detail  through  the  use  of  dis¬ 
tribution  fvinctions.  It  is  shown  that  the  functions  are  discontinuous  in 
velocity  variables  and  lead  to  the  same  results  as  particle  trajectory 
theory.  The  limitations  of  the  sheath  theory  and  its  connection  with 
space-charge -limited  diode  theory  are  discussed.  It  is  shown  that,  under 
certain  assumptions,  the  potential  on  a  satellite  whose  diameter  is  smaller 
than  the  local  Debye  length  will  reach  3500  volts  negative  in  the  more 
Intense  regions  of  the  Van  Allen  electron  belt.  The  equilibrixim  potential 
decreases  with  increasing  size  of  the  satellite,  with  a  limiting  value  of 
-35,000  volts  for  satellites  much  greater  than  the  Debye  length  in  diameter. 
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I.  INTRODUCTION 


The  initial  flow  of  electrons  to  an  absorbing  body  placed  in  a  neu¬ 
tral  isothermal  plasma  will  exceed  the  flow  of  positive  ions,  and  the  body 
will  begin  to  accumulate  a  net  negative  charge.  As  the  charge  builds  up, 
It  will  repel  more  plasma  electrons,  attract  more  ions,  and  the  charging 
rate  will  tend  toward  zero. 

A  space  vehicle  moving  through  the  plasma  generally  assumed  to  exist 
at  and  beyond  several  earth  radii  from  the  earth  should  tend  to  become 
negatively  charged  in  a  similar  manner,  although  the  amount  of  charge  ac¬ 
cumulated  may  be  considerably  modified  by  solar  photoeffects  and  secondary 
emission.  The  intense  flux  of  geomagnetically  trapped  electrons  in  the 
outer  radiation  belt  should  produce  an  unusually  large  negative  charge  on 
a  vehicle  penetrating  this  region. 


Solar  coronal  theories  indicate  that  the  interplanetary  gas  in  the 
region  of  the  earth's  orbit  consists  primarily  of  protons  and  electrons 
at  temperatures  of  at  least  20,000°  K,  and  that  the  coronal  gas  streams 
out  from  the  sun  with  its  stream  velocity  large  compared  with  the  proton 
mean  thermal  speed.  This  solar  wind  can  normally  penetrate  the  sunward 


side  of  the  earth's  magnetic  field  to  a  distance  of  about  fifteen  earth 
radii  from  the  earth's  center.  The  condition  of  the  gas  within  the 
shielded  region  (the  magnetosphere)  around  the  earth  is  not  well  known. 
Chapman, for  example,  suggests  the  temperature  above  the  earth's  atmos¬ 
phere  rises  gradually  to  the  outside  value,  whereas  Johnson  proposes  the 
temperature  throughout  the  magnetosphere  is  relatively  constant  at  about 
1250  K.^  ^  Whistler  data  are  reasonably  consistent  with  the  higher- tem¬ 

perature  model,  while  satellite  ion  probe  data  tend  to  show  much  lower 
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temperatures.  Neither  model  can  claim  confirmation  at  this  time. 

Satellite  charge-up  in  the  ionosphere  has  been  treated  by  Jastrow  and 
Pearse,'  who  assume  the  satellite  potential  <f)^  is  negative  and  that  the 
satellite  is  surrounded  by  a  thin,  positively  charged  sheath.  With  the 
satellite  velocity  intermediate  between  the  ion  and  electron  velocities, 
they  find  the  ion  density  around  the  satellite  unchanged  from  its  ambient 
ionospheric  density,  whereas  the  electron  density  is  decreased  by  the  elec¬ 
tron  Boltzmann  factor  e"^,  where 


Under  these  conditions,  the  ion  flux  at  the  satellite  is  equal  to  the  pro¬ 
duct  of  the  ambient  ion  density  and  the  satellite  velocity  v^,  and  the 
electron  flux  is  equal  to  the  product  of  the  ambient  electron  flux  and  the 
Boltzmann  factor.  The  ambient  flux  is  one-fourth  the  product  of  the  am¬ 
bient  density  and  mean  thermal  velocity  v^,  and,  when  the  electron  flux 
equals  the  ion  flux. 


X  = 


In  their  .jag  calculations,  Jastrow  and  Pearse  find  that  the  effective 
cross  section  for  ion- satellite  collisions  is  about  twice  the  satellite 
geometric  cross  section.  Their  neglect  of  this  effect  (which  will  decrease 
X  by  In  2)  in  the  potential  calculation,  and  their  somewhat  incorrect  form 
for  the  electron  flux  (their  Eq.  (4)  results  in  a  slightly  modified  Boltz¬ 
mann  factor),  will  not  alter  significantly  their  resrrlt  that  xss  7. 
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Tr»  contrast,  Chopra^ conjectures  that  the  shielding  of  the  Coulomb 
potential  of  the  satellite  is  negligible,  from  which  he  concludes  that  the 
potential  of  a  body  at  rest  in  an  ionized  gas  is  given  by 


V 


where  is  the  ion  meein  thermal  velocity.  It  can  be  shown  (see  Section 
VI  following)  that  the  above  result  for  x  should  be  decreased  by  ln(l  +  x). 
It  is  shown  in  Section  VI  that  Chopra's  result  is  correct  for  conditions 
opposite  those  he  asstimed,  and  is  applicable  to  a  body  completely  shielded 
by  a  vanishingly  thin  sheath.  For  a  plasma  of  electrons  and  protons  at 
equal  temperat\ires,  it  is  shown  that  the  extreme  values  are  x  =  2.5 
shielding)  and  x  =  3.8  (con5)lete  shielding) , and  that  intermediate  shielding 
gives  Intermediate  values  of  x. 

Although  details  of  the  sheath  structure  are  apparently  unnecessary 
for  csCLculatlon  of  the  potential  on  the  satellite,  they  have  impoirtant 
bearing  on  drag  and  other  effects.  If  the  sheath  is  sufficiently  thin  and 
the  satellite  sufficiently  attractive  to  the  ions,  it  is  sometimes  pos¬ 
sible^  to  make  order-of-magnitude  estimates  of  the  charge  density  in  the 
sheath  with  the  aid  of  conventional  probe  theory.  Probe  theory  in  the 
usual  sense,  however,  often  coiiq>ounds  the  inconsistencies  of  the  idealized 
sheath  with  generally  inapplicable  diode  theory,  and  is  of  little  general 
value.  More  recent  probe  theories,  such  as  that  of  Bernstein  and  Rabino- 
vitz,^^^  attempt  to  determine  the  particle  distribution  functions  from 
Boltzmann's  equations.  Such  treatments  encounter  many  mathematiceO.  diffi¬ 
culties,  however,  and  no  general  solutions  have  been  obtained.  Ichimlya, 
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(7) 

et  al.,  in  their  discussion  of  rocket  ion  probes,  assume  (see  their 
Eqs.  (1)  and  (2))  that  the  ion  current  to  the  probe  is  given  both  by  the 

/  Q\ 

Mott-Smith  expression  based  on  the  ideal  sheath  condition,  and  by  the 
expression  for  the  current  in  a  cold  spherical  diode.  Although  the  latter 
expression  can  be  shown  to  be  the  zero  temperature  limit  of  the  former 
(under  certain  further  restrictions;  see  the  discussion  at  the  end  of  Sec¬ 
tion  IV),  they  equate  the  two  expressions  in  order  to  eliminate  the  sheath 
thickness.  Similarly,  Hundley^^^  relies  heavily  on  diode  theory  to  calcu¬ 
late  the  satellite  potential  in  the  Van  Allen  regions,  and  finds  it  neces¬ 
sary  to  attempt  to  approximate  the  cumbersome  diode  solutions  with  simpler 
e3q)resslons.  However,  it  will  be  shown  that  even  in  the  presence  of  the 
Van  Allen  electrons,  potentials  calculated  on  the  basis  of  the  ideal  sheath 
theory  vary  only  by  a  factor  of  two  or  three  from  an  intermediate  value 
over  the  entire  range  of  possible  sheath  thicknesses.  Since  the  errors 
from  the  uncertainties  in  the  characteristics  of  the  ambient  plasma  and 
Van  Allen  flux  cause  considerably  more  variation,  Hundley*  s  sheath  thick¬ 
ness  estimates,  however  incorrect,  are  largely  irrelevant  to  his  estimate 
of  the  satellite  potential. 

It  is  the  purpose  of  the  following  to  demonstrate  the  simplicity  of 
the  ideal  sheath  theory  in  the  evaluation  of  the  current  to  a  body  at  rest 
in  collisionless  ionized  gas,  and  to  apply  these  results  to  the  satellite 
potential  in  the  Van  Allen  region.  In  view  of  the  occasional  confusion  on 
the  subject,  and  in  order  to  provide  a  general  foundation  for  the  theory, 
the  method  of  the  distribution  function  and  the  Boltzmann  equation  is  used. 
It  is  shown  that  certain  discontinuous  distribution  functions  lead  to  the 
same  result  as  the  particle  trajectory  methods,  so  that,  in  pai-ticvaar, 
the  results  of  Mott-Smith  and  Langmuir^®)  are  recovered  and  extended. 
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II.  METHOD 

In  order  to  estimate  the  potential  on  a  satellite  in  the  radiation 
belt,  several  simplifying  assumptions  will  be  made.  It  will  be  assumed 
that  the  gas  is  an  infinite,  homogeneous,  collisionless  plasma  of  protons 
and  electrons  with  isotropic  velocity  distribution,  and  that  nearly  all 
proton  velocities  are  much  greater  than  the  satellite  velocity.  Since  the 
velocity  of  a  satellite  in  the  outer  radiation  belt  (at  about  four  earth 
radii  from  the  center  of  the  earth)  equals  the  average  proton  speed  in  a 
Maxwellian  distribution  at  1000°  K,  this  assumption  is  very  good  for  the 
continuously  rising  temperature  assumption,  but  poor  for  the  constant  low 
temperature  assumption.  The  high  temperature  approximation  assures  spher¬ 
ical  symmetry  to  the  problem  and  permits  analytic  evaluations  of  the  pro¬ 
ton  currents  to  the  satellite.  Furthermore,  the  results  so  obtained  will 
be  directly  applicable  to  static  floating  probe  theory.  If  the  satellite 
and  average  proton  speeds  are  assumed  comparable,  the  symmetry  is  lost 
and  numerical  methods  will  be  required.  Even  if  the  proton  temperature  is 
as  low  as  1250  K  the  results  based  on  the  high  temperature  approximation 
will  provide  reasonable  estimates  of  the  satellite  potential,  although 
many  of  the  details  of  the  proton  current  calculations  will  be  more  pre¬ 
cise  than  the  approximation  warrants. 

The  perturbing  effects  of  solar  photons,  etc.,  will  be  neglected^  the 
satellite  will  be  assumed  to  be  a  conducting  sphere  which  absorbs  all  inci¬ 
dent  particles,  and  the  electrons  will  be  treated  non-relatlvlstically. 
Although  the  only  quantity  to  be  calculated  is  the  potential  of  the  satel¬ 
lite,  it  is  of  interest  to  introduce  the  proton  and  electron  distribution 
functions  and  the  equations  they  satisfy,  the  collisionless  (homogeneous) 
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Boltzioann  equations.  If  the  origin  of  the  coordinate  system  is  at  the 
center  of  the  satellite,  then,  by  symmetry,  the  electron  ajid  proton  cur¬ 
rents  must  be  radial,  and  each  distribution  function  may  be  divided  into 
two  parts  corresponding  to  inward — or  outward- -moving  particles.  Usually, 
in  the  steady  state,  only  the  inward  moving  particles  will  be  of  interest, 
and  it  is  convenient  to  identify  the  flux  of  such  particles.  The  expression, 
"inward  electron  flux,"  for  example,  will  refer  to  the  flux  of  electrons 
with  negative  radial  velocity  components.  The  radial  electron  current  is 
the  difference  between  the  outward  and  inward  electron  flux.  In  an  infi¬ 
nite  uniform  gas  with  an  isotropic  velocity  distribution,  the  flux  in  any 
direction  is  one- fourth  the  product  of  the  gas  density  and  the  mean  par¬ 
ticle  speed,  and  the  current  is  zero. 

One  boixndary  condition  on  the  solutions  to  the  Boltzmann  equation, 
that  of  the  perfectly  absorbing  sphere,  means  that  the  outward  electron 
flux  and  outward  proton  flux  each  vanish  at  the  surface  of  the  sphere. 

Another  boundary  condition  characterizes  the  plasma  at  infinity  or,  more 
generally,  beyond  a  particular  radial  distance  (possibly  infinite)  from 
the  center  of  the  satellite.  It  will  be  assumed  that  at  this  radial,  dis¬ 
tance,  which  will  be  called  the  sheath  edge,  the  inward  proton  and  inward 
electron  distribution  functions  are  known,  and  are  isotropic  in  their  ve¬ 
locity  variables.  Essentially,  the  isotropic  assumption  implies  the  usual 
idealized  sheath  conditions,  i.e.,  that  the  potential  outside  the  sheath 
region  vanishes,  and  that  the  plasma  there  is  in  an  unperturbed  state.  The 
assumption  is  therefore  contradictory  (since  in  the  presence  of  the  satel¬ 
lite  there  will  be  radial  currents  of  electrons  and  protons  throughout  the 
entire  plasma)  for  a  finite  sheath,  although  hopefully  it  is  a  good 
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approximation  if  the  potential  at  the  edge  of  the  sheath  is  sxifficiently 
small. 

Since  the  distribution  functions  vanish  inside  the  absorbing  sphere, 
they  satisfy  the  Boltzmann  equations  trivially  in  this  region.  The  sur¬ 
face  of  the  sphere  represents  a  surface  of  discontinuity  in  configuration 
space  for  the  distribution  functions,  and  a  surface  on  which  the  Boltzmann 
equation  does  not  apply.  Similarly,  in  the  presence  of  the  negative  poten¬ 
tial,  there  are  regions  in  phase  space  where  the  sum  of  the  proton  kinetic 
and  potential  energies  is  negative,  and  throughout  which  the  proton  distri¬ 
bution  function  vanishes  in  the  absence  of  collisions.  The  surface  of  zero 
proton  total  energy  is  a  surface  of  discontinuity  for  the  proton  distribu¬ 
tion  function,  and  a  surface  on  which  the  Boltzmann  equation  is  not  satisfied. 

Since  the  incoming  (outgoing)  electron  distribution  function  vanishes 
for  positive  (negative)  radial  velocities,  the  zero  radial  velocity  sur¬ 
face  in  electron  velocity  space  is,  in  general,  a  surface  of  discontinuity 
for  the  partial  electron  distribution  functions.  In  the  steady  state,  the 
electron  partial  distribution  functions  independently  satisfy  the  Boltz¬ 
mann  equation  everywhere  except  on  this  surface.  Similar  considerations 
hold  for  the  proton  partial  distribution  functions,  except  that  the  sur¬ 
face  of  discontinuity  is  the  zero  totail  energy  surface  plus  that  section 
of  the  zero  radial  velocity  surface  which  lies  in  the  positive  total  energy 
region  of  proton  phase  space. 

Steady  state  solutions  for  proton  and  electron  inward  moving  distri¬ 
bution  functions  will  be  obtained.  Each  inward  moving  distribution  func¬ 
tion  vanishes  on  one  side  of  its  surface  of  discontinuity;  on  the  other  side 
it  is  a  solution  of  the  collisionless  Boltzmann  equation  which  satisfies 
the  boundary  condition  for  the  region  outside  the  sheath. 
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III.  THE  ELECTRON  FLUX 

The  ccllisionless  Boltzmann  equation  states  that  the  distribution 
function  is  constant  along  a  particle  trajectory.  Assume  the  electrostatic 
potential  0^  at  the  surface  of  the  satellite  is  negative  and  that  0(r)  >  0 
for  all  r  >  R,  where  R  is  the  radius  of  the  satellite.  The  velocity  dis¬ 
tribution  of  inward  moving  particles  is  assumed  isotropic  for  r  >  a,  where 
a  is  the  sheath  radius.  For  r  >  a  the  electron  partial  distribution  func¬ 
tion  f  IS  a  function  only  of  the  kinetic  energy.  Since  the  total  ener¬ 
gy  of  an  electron  is  conserved  along  its  trajectory,  the  (inward)  distri¬ 
bution  function  Is  f(m  +  z),  where  z  =  -  e0^  and  is  the  potential  energy 
of  an  electron  at  R. 

The  flux  1  at  a  point  r  in  a  given  direction  may  be  caJ.culated  from 
the  partial  distribution  function  for  particles  moving  in  that  direction. 
With  spherical  coordinates,  the  velocity  component  in  the  given  direction 
IS  V  cos  «,  and  the  flux  will  involve  an  integration  over  velocity  space 
of  this  factor  times  the  partial  distribution  function.  When  the  velocity 

distribution  is  isotropic,  the  integration  over  angles  may  be  carried  out 
to  give 


r  = 


r)  do) 


(1) 


where  A  is  a  constant.  In  the  case  of  the  radially  inward  moving  electrons 
in  the  spherically  symmetric  field,  the  inward  flux  at  R  is  given  by 


Jo 


ajf(<jo  +  z)  do) 


(2) 


The  inward  fliix  for  r  >  a  is  given  by 


r 

o 


V 

e 


(3) 


where  is  the  flux  in  any  direction  in  the  absence  of  the  satellite,  N 
is  the  density,  and  the  mean  speed  of  the  unperturbed  electrons. 

If  the  partial  distribution  function  for  r  >  a  is  Maxwellian,  that  is, 

if  it  goes  as  exp(-a3/kT),  then  for  r  =  R,  it  is  changed  by  the  factor  e"^ 
where 


X 


kT  kT 


(^) 


The  factor  e  may  be  brought  outside  the  integral  sign  in  Eq.  (2),  and  it 
follows  immediately  that  the  inward  electron  flux  to  the  satellite  is 


Instead  of  the  distribution  function,  the  "integral  flux"  t  (E)  above 
the  energy  E  may  be  given  as  a  boimdary  condition,  where 


i(E)  =  A 


CO 

a)f(cD)dto 


(6) 


and  \|r(0)  =  From  Eq.  (6), 
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where  Lhe  prime  indicates  differentiation.  Substitution  of  Eq.  (7)  into 
Eq.  (2)  gives  the  inward  electron  flux  to  the  satellite  in  terms  of  the  un¬ 
perturbed  integral  flux  as 


r 


03 


=  'l'(z)  + 


doj 


Integration  by  parts  gives 


du> 


(8) 


CO 


An  example  of  the  integrea  flux  is  the  monoenergetic  source  at  energy 

E  ,  that  is,  where 
o' 


E  <  E 
—  o 


t(E)  =<  ° 

0  E  >  E 


where  T  is  a  constant.  From  Eq.  (8)  the  flux  at  the  satellite  becomes 
o 


(1  -  f)  fo 

o 


z  <  E 
—  o 


z  >  E 


0 
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For  sufficiently  high  energies,  the  integral  flux  in  the  outer  belt 
of  trapped  electrons  is  assumed  to  decrease  inversely  as  some  power  of  the 
energy.  If,  for  energies  E  >  z,  the  unperturbed  flux  i^(E)  goes  as  E 
with  n  >  0,  then,  from  Eq.  (8),  the  flux  at  the  satellite  becomes 

n  +  1 


(9) 
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IV.  THE  PROTON  Pl-UX 

The  (inwai'd)  proton  distribution  function  may  be  treated  in  the  same 
manner  as  vers  the  electron  functions  except  for  those  regions  of  phase 
space  where  the  sum  of  the  kinetic  energy  u)  and  the  potential  energy  e0j^  =  -z 
is  less  than  zero.  In  these  "forbidden"  regions  where  cu  <  z  the  distribu¬ 
tion  function  vanishes,  since  an  originally  free  particle  cannot  (in  the 
absence  of  collisions)  move  into  a  region  where  its  total  energy  is  nega¬ 
tive. 

In  order  to  find  the  forbidden  region  of  velocity  space  at  the  sur¬ 
face  of  the  satellite,  consider  a  proton  with  radia.1  velocity  u  (positive 

inward)  and  tangential  velocity  t  at  the  surface.  If  u  and  t  are  the 

o  o 

corresponding  velocities  at  the  edge  of  the  sheath,  then,  by  conservation 
of  energy. 


2 

u  + 


+ 


M 


(10) 


where  M  is  the  proton  mass.  (The  relation  in  Eq.  (lO)  is  a  necessary  con¬ 
dition  for  a  particle  reaching  the  satellite.  The  sufficient  condition, 
that  is,  the  condition  on  the  potential  in  the  intermediate  region  such 
that  all  particles  satisfying  the  conservation  laws  actually  reach  the  sur¬ 
face  will  be  discussed  later.)  Conservation  of  angular  momentum  gives 

Rt  =*at^  (11) 


and  may  'be  eliminated  from  Eqs.  (10)  and  (ll)  to  give 


di 

u  + 


M 


(12) 
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If,  as  in  a  Maxwellian  plasma,  a  continuous  range  of  initial  radial  speeds 
from  zero  to  infinity  exists,  the  distribution  function  at  R  vanishes  iden¬ 
tically  if,  and  only  if,  the  corresponding  velocity  variables  are  such  that 
the  expression  given  by  the  left  side  of  Eq.  (12)  is  negative. 

A  convenient  way  to  calculate  the  proton  flux  at  the  satellite  surface 
is  to  assume  at  first  that  if  the  incoming  proton  distribution  function  is 
f  (oj)  at  r  =  a,  then  it  is  given  by  f  (tn  -  z)  at  r  =  R  for  positive  o), 

and  to  subtract  the  flux  associated  with  regions  where  a>  <  z.  For  a  Max¬ 
wellian  distribution,  it  follows  immediately  that  the  flux  of  protons  is 
equal  to  e^  times  the  difference  between  the  unperturbed  values  of  the 
total  and  the  "forbidden”  flux,  that  is. 


r 

P 


(13) 


where  the  subscript  p  refers  to  protons.  The  ratio  of  the  forbidden  to 
total  flux  is 


where,  from  Eq.  (12),  the  upper  limits  for  the  forbidden  region  are  given 

by 


2 


U 


) 


2  _  _  2z 

^1  ■  ■  - W 


(14) 
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The  integration  may  be  carried  out  and  the  results  expressed  as 


r  =  T- 
P  po 


(15) 


where 


(16) 


is  a  positive  parameter  characterizing  the  relative  sizes  of  the  satellite 
and  sheath.  This  resvLLt  was  first  obtained  on  the  basis  of  particle  tra¬ 
jectory  methods  by  Mott-Smith  and  Langmuir.  A  satellite  potentiad  may 
be  regarded  as  strong  if  »  1,  that  is,  if 


X  =  -  -  1 

kT  >>  ^2  -  n 


(17) 


for  then  the  proton  flux  at  the  satellite  approaches 


r  — r  I— — ..ii  1  =  r  — 

P  PO\  u  /  po  p2 


(18) 


This  limiting  value  may  also  be  interpreted  as  the  thin  sheath  approxima- 
tion,  eilthough  strictly  it  is  always  necessetry  to  consider  both  potential 
and  sheath  size  together.  In  the  opposite  limit,  ux  «  1,  that  of  weak 
potentials  or  thick  sheath,  the  flux  approaches 


r 

po 


(1  +  x) 


r 

po 


kT  -  e0j^ 
kT 


(19) 
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V.  COMSTRAIMTS  ON  POTENTIAL  SHAPES 

The  preceeding  results  are  valid  only  if  the  potential  throTighout  the 
sheath  satisfies  certain  conditions,  and  under  more  general  conditions  the 
calculations  may  proceed  stepwise.  The  conditions  are  the  same  for  all 
methods  based  on  the  use  of  energy  and  angular  momentum  conservation  at  two 
separated  points,  and  are  imposed  to  insure  that  the  proper  trajectory  con¬ 
necting  the  two  points  exists. 

The  condition  on  the  validity  of  Eq.  (5)  for  electrons  moving  in  a 
repulsive  field  is  the  rather  weak  restriction  that  throughout  the  sheath 
region 

0(r)  >  (20) 


This  condition  is  easily  deidved  when  the  time- reversed  problem  of  an  elec¬ 
tron  directed  from  the  surface  of  the  satellite  into  the  sheath  is  consid¬ 
ered.  If  electron  speeds  down  to  zero  are  permitted,  then  if  and  only  if 
Eq.  (20)  holds  will  all  ejected  electrons  escape  the  sheath. 

The  condition  for  the  validity  of  Eq.  (15)  for  the  proton  flux  may  be 
obtained  from  the  requirement  that  protons  which  penetrate  the  entire  sheath 
do  not  reach  zero  radial  velocity  in  the  sheath.  From  the  conservation  laws 
this  requirement  will  be  met  thro\ighout  the  sheath  only  if 


2  ^2.  Ma  .Z 


(21) 


where  0  is  the  potential  at  r,  and  u^  and  t^  are  the  velocity  conrponents  at 
a,  the  sheath  edge.  If  u,  and  t  are  the  corresponding  values  at  R,  the 
surface  of  the  satellite,  then 
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M 

2 


2 


u  + 


(22) 


2 

Equation  (22)  may  be  solved  for  and  the  result  substituted  into  Eq.  (21) 
to  give 


(23) 


As  u  approaches  zero  and  t  approaches  t^,  given  by  Eq.  (l4),  u^  also  ap¬ 
proaches  zero,  and  the  two  positive  terras  on  tne  rignt-nand  side  of  Eq.  (23) 
can  vanish  simultaneously.  It  is  therefore  necessary  and  sufficient  that 
for  all  R  <  r  <  a, 

0R  (24) 


0(r)  <  0j^(r) 


2 

r 


-  1 


in  order  that  the  proton  flux  be  given  by  Eq.  (15).  The  limiting  potential 

0  (r)  is  proportional  to  l/r  plus  a  constant,  with  the  constants  deter- 
M 

mined  by  0w(r)  =  0„  and  0j,(a)  =  0.  The  restriction  of  the  potential  to 
R  M, 

values  less  than  fi),,  will  be  called  the  sheath  condition.  In  the  absence 
M 
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of  infinite  charge  densities,  solutions  of  Poisson's  equation  which  vanish 
for  r  >  a  cannot  satisfy  the  sheath  condition,  but  if  the  point  where  the 
potential  rises  above  0^^  is  near  the  sheath  edge  the  error  in  Eq.  (15)  will 
be  correspondingly  small. 

If  the  results  in  Eqs.  (5)  and  (15)  are  simultaneously  valid,  the  po¬ 
tential  must  satisfy 

^  <  0  <  0„ 

A  more  restrictive  lower  limit  may  be  obtained  from  the  reasonable  assump¬ 
tion  that  the  charge  density  is  positive  throxighout  the  sheath.  If  this 
condition  holds,  the  electric  field  intensity  inside  the  sheath  is  nega¬ 
tive, 

0'  >  0 

and  0  Increases  monotonically  thro\ighout  the  sheath.  If  0  is  a  potentieG. 

which  varies  as  l/r  plus  a  constant  and  satisfies  0„(R)  =  0^  and  0_(a)  =  0, 

C  R  C 

then 


0c"  -  -  I  0c'  (») 

that  is,  0  represents  a  "sheath"  of  zero  charge  density.  In  a  positively 
charged  sheath  it  is  necessary  that  0  satisfy  the  inequEility 

0"  <  -  I  0. 

Since  0  is  a  solution  of  Poisson's  equation,  0  is  greater  than  0  near  the 
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sheath  edge,  and  it  can  be  shovm  from  Eq,  (25)  that  the  cuirves  represent¬ 
ing  0  and  0^  cannot  cross  Inside  the  sheath.  For  if,  as  r  decreases  from 
the  sheath  edge,  the  potential  0  falls  below  0  ,  its  slope  must  continue 
to  decrease  so  rapidly  in  order  to  satisfy  Eq.  (25)  that  0(r)  will  be  be- 
low  0j^.  Hence,  everywhere  inside  the  positively  charged  sheath,  the  po¬ 
tential  0  must  be  greater  than  0^.  This  lower  limit  on  the  potential  is 
more  restrictive  and  can  replace  that  given  by  Eq.  (20),  and  lnq)lles  that, 
if  the  proton  flux  calciilation  is  valid  in  a  positively  charged  sheath,  the 
potential  inside  the  sheath  is  restricted  to  the  rather  narrow  range 

0c  <  l<  <  i#„ 
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VI.  THE  SATEILLITE  POTENTIAL 

In  a  neutral  Maxwellian  plasma  the  satellite  will  attain  a  potential 
such  that  the  electron  flux  in  Eq.  (5)  equals  the  proton  flux  in  Eq.  (15) 
at  the  siirface,  that  is. 


r  e-^  =  r 


po 


1  + 


(26) 


The  electron  and  proton  temperatures  have  been  ass\amed  equal  in  Eq.  (26), 
although  a  ten^erature  ratio  factor  may  be  introduced  in  one  of  the  ex¬ 
ponents  if  necessary.  With  equal  teii5)eratures  and  number  densities,  the 
unperturbed  electron- to-proton  flux  ratio  is  inversely  proportional  to  the 
square  root  of  the  corresponding  mass  ratio,  and  Eq.  (26)  becomes 


1/2 


-px  .  /M  '  -X  , 
e  +  p  -j  e  =  1  + 


(27) 


If  the  sheath  thickness  is  zero  (a  =  R,  ^  Eq.  (2?)  yields 


or  X  =  3.8  for  protons  and  electrons.  If  the  sheath  is  infinite  (a  =  », 
p  =  0),  Eq.  (27)  becomes 


X  .  i„(a  .  X)  .  i  1„  (I) 

and  X  =  2.5.  From  Eq.  (27)  the  satellite  potential  increases  as  a  in 
creases  and  0^  lies  between  -3.8  kT/e  and  -2.5  kT/e. 
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If  the  electron  flux  is  expressed  in  terms  of  the  integral  flux  i|f(E) 
defined  by  Eq.  (6),  then,  from  Eq.  (8),  the  equilibrium  condition  is 


00 

r  . 

J  dE  =  r 

1  +  u  -  e 

-4  po 

u 

(28) 


where  z  -  xkT.  For  example,  the  integral  flux  for  the  plasma  electrons  is 


i(E) 


(29) 


so  that  when  Eq.  (29)  is  substituted  into  Eq.  (28)  and  the  integration  per¬ 
formed,  Eq.  (26)  is  recovered.  In  the  region  of  the  outer  radiation  belt, 
the  flux  of  plasma  electrons  to  the  charged  satellite  is  negligible  in  com¬ 
parison  with  the  flux  of  high-energy  trapped  electrons,  and  the  expression 
in  Eq.  (28)  is  more  convenient. 

If,  for  energies  E  >  z,  ^  goes  as  E""^,  then,  from  Eq.  (9)  Eq.  (28)  be- 

comes 


i(z) 

n  +  1 


r 

po 


(30) 


If  z^  and  z^  represent  the  values  of  z  for  zero  and  infinite  sheath  thick¬ 
ness,  then,  from  Eq.  (30) 


\|f(z  )  =  (n  +  1)  r 
o  po 


(31) 
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and 


00 


(n  -t-  i) 


po 


(32) 


since  a  >>  kT. 

If  the  proton  density  in  this  region  is  lOO/cm^  and  the  plasma  tem¬ 
perature  is  5000°  K,  the  unperturbed  proton  flux  f  is  approximately 
7  2 

10  protons/cm  /sec.  Since  n  is  usually  considered  to  be  in  the  range  3  to 

^  the  value  of  is,  from  Eq.  (31) ^  that  energy  above  which  the  un- 

7  2 

perturbed  electron  flux  is  ~  5  x  10 '  electrons/cm  /sec.  In  the  heart  of 
the  outer  belt,  this  flux  occurs  above  an  energy  of  about  55  kev,^^^^  al¬ 
though  occasionally  a  larger  flux  is  reported. ^ An  increase  in  flux  by 
a  factor  of  100  (with  the  same  power  law)  increases  by  a  factor  of  three. 

The  satellite  potential  increases  as  the  sheath  thickness  increases. 
From  Eqs.  (3I)  and  (32),  z  differs  from  z  by  the  factor 

00  O 


(33) 


Thus,  the  potential  for  the  infinite  sheath  is  about  l/lO  that  for  zero 
sheath,  and  again,  the  potential  is  not  strongly  dependent  on  the  thickness. 

As  was  pointed  out  in  Eq.  (17)  it  is  not  sufficient  to  consider  only 
the  relative  size  of  the  sheath  and  satellite  in  order  to  estimate  the  po¬ 
tential  as  that  given  by  the  thick  or  thin  sheath  limit.  The  relevant  cri¬ 
terion  is  whether  |ix  is  much  larger  or  much  smaller  than  unity.  Since 
X  »  1,  it  follows  from  Eq.  (I7)  that  a  sheath  is  thin  (thick)  when  its 
radius  is  small  (large)  compared  with  Rx^/^.  If  x  =  10^  for  thick  and 
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5 

X  =  10  for  thin  sheaths,  then  z  is  a  good  approximation  when  a  >  >  100  R, 

v;hile  z  holds  rflien  a  <  <  jiOO  R. 
o 

In  order  to  determine  the  sheath  thickness  and  test  the  consistency  of 
the  model,  it  is  necessary  to  determine  the  outgoing  distribution  functions 
over  velocity  space  in  order  to  obtain  an  expression  for  the  charge  density 
as  a  function  of  the  potential.  Poisson's  equation  may  be  used  to  elimi¬ 
nate  the  charge  density,  and  the  resulting  equation  for  the  potential  must 
be  solved  subject  to  the  necessary  boundary  conditions.  In  general,  this 
becomes  very  involved  mathematically,  and  only  approximate  solutions  to 
certain  special  cases  have  been  obtained. 

Accurate  solutions  may  be  obtained  in  the  limiting  case  where  the  elec¬ 
tron  and  proton  energies  in  the  undisturbed  plasma  approach  zero  while  their 
densities  approach  infinity  in  such  a  manner  that  the  undisturbed  flux  of 
each  is  finite.  In  this  case  the  ingoing  and  outgoing  electron  distribu¬ 
tion  functions  and  the  outgoing  proton  distribution  function  vanish  in  the 

sheath.  The  proton  density  will  be  fovind  to  be  inversely  proportional  to 
2  l/  2 

r  (-0)  f  and  a  solution  (in  the  form  of  an  infinite  series)  to  Poisson's 
equation  may  be  found  for  which  the  potential  and  its  first  derivative 
vanish  at  the  sheath  edge.  These  solutions  were  first  derived  in  the  study 
of  space- charge-limited  diodes  with  cold  emitters  where  the  charge  density 
between  the  diode  elements  follows  directly  from  current  and  energy  conser¬ 
vation  laws. 

It  might  seem  possible  to  make  order  of  magnitude  estimates  of  the 
sheath  thickness  aivund  a  satellite,  based  on  the  cold  emitting  space- 
charge-  limited  diode  if  the  emitting  surface  is  imagined  inside  the  sheath 
at  a  radius  where  the  potential  Is  sufficiently  strong  so  that  very  few 
electrons  are  present  and  the  protons  are  moving  Inward  primarily  in  the 
radial  direction.  This  approach  has  had  some  success  in  correlating  gas 
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discharge  measurements  on  Maxwellian  plasmas  with  plane  probes.  Since  in 
this  case  the  electron  density  probably  decreases  exponentially  with  the 
Boltzmann  factor,  it  may  be  neglected  wherever  -e0  is  a  few  times  kT,  while 
energy  considerations  show  that  the  proton  velocities  in  these  regions  tend 
to  be  largely  radial  in  direction.  However,  the  emitter  is  no  longer  cold, 
since  the  protons  have  nonzero  initial  velocities;  the  emitter  is  not  space- 
charge-limited,  since  an  (unknown)  electric  field  exists  at  the  emitter; 
and  the  emitter  current  is  not  known.  In  the  case  of  the  satellite,  the 
electron  density  in  the  sheath  does  not  decrease  as  rapidly  with  potential, 
because  of  the  presence  of  the  high-energy  electrons.  It  has  been  assumed, 
however,  that  the  density  of  high-energy  electrons  in  the  radiation  belt 
was  at  least  an  order  of  magnitude  less  than  that  of  the  Maxwellian  elec¬ 
trons,  and  therefore  their  density  in  the  sheath  may  be  neglected,  (it 
should  be  noted  that  this  assumption  requires  that  the  steep  power  law  de¬ 
pendence  of  the  high-energy  flux  cannot  hold  for  energies  much  below  10 
kev. ) 

In  order  to  avoid  a  net  negative  charge  density  in  the  sheath,  the 
potential  cannot  decrease  too  rapidly  inward  from  the  sheath  edge,  and 
the  imagined  emitter  surface  may  be  well  inside  the  sheath.  At  best,  it 
might  be  hoped  that  sheath  thickness  calculations  based  on  space- charge- 
limited  diode  theory  give  an  order  of  magnitude  estimate  of  the  lower  limit 
of  the  sheath.  With  these  reservations,  diode  theory  may  be  shown  to  be 
consistent  with  the  estimate  that  a  sheath  is  large  (small)  conrpared  with 

1/2 

Rx  when  the  satellite  radius  R  is  small  (large)  conq)ared  with  the  Debye 
length  D  of  the  undisturbed  plasma  (D  =  40  cm  for  the  density  and 
teiiq)erature  given) . 
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VII.  IMPACT  PARAMETER  METHODS 

Although  extensions  of  the  preceding  calculations  are  most  convenient¬ 
ly  carried  out  through  the  use  of  the  distribution  functions,  these  pre¬ 
liminary  results  can  be  obtained  in  a  somewhat  more  physical  method  through 
considerations  of  the  impact  parameters  for  grazing  collision.  It  follows 
immediately  from  conservation  of  angular  momentum  and  energy  that  all  par¬ 
ticles  with  impact  parameters  less  than 


R 


'E  -  q^. 
.  E 


where  E  is  the  particle  energy  outside  the  sheath  where  0=0,  and  q  is 
the  particle  charge,  will  strike  the  satellite  of  radius  R  and  potential 
0p,  and  the  satellite  cross  section  for  collision  with  a  particle  of  charge 
q  is  equal  to  its  geometrical  cross  section  times  the  factor 


^  -  ‘^R 
E 


These  relations  hold,  only  if  the  sheath  condition  in  Eq.  (24)  holds,  and, 
in  addition,  the  impact  parameters  lie  between  zero  and  a. 

If  the  Integral  electron  flvix  outside  the  sheath  is  i|r(E),  the  fltix  of 
electrons  with  energies  between  E  and  E  +  dE  is 

-d\|f  =  -i|r*(E)  dE 


and  the  flux  at  the  satellite  becomes 
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r  = 


\|f'(E)  dE 


(37) 


where  q  =  -e.  The  lower  limit  is  -e0„,  since  a  particle  with  initial  ener- 

R 

gy  E  <  -e0  cannot  penetrate  to  the  satellite. 

R 

Integration  by  parts  gives 


r 


dE 


which  is  identical  to  Eq.  (8),  the  result  obtained  with  the  use  of  the  dis¬ 
tribution  function. 

If  the  density  in  velocity  space  is  given  instead  of  the  integral  flux, 
the  differential  flux  may  be  written  directly  in  terms  of  the  density.  For 
an  isotropic  distribution,  the  number  of  particles  with  speeds  between  v 

p 

and  V  +  dv  is  proportional  to  n(E)v  dv,  where  n(E)  is  the  velocity  space 

2 

density.  The  differential  flux  is  then  vn(E)v  dv,  or,  with  the  proper 
normalizing  factor  incorporated  into  n(E), 


-di/  =  En(E)  dE 


(38) 


This  expression  may  be  substituted  into  Eq.  (37)  give 


r  = 


(E  +  e0  )  n(E)  dE 
K 


or 


r  =y  En(E  -  dE 


(39) 


(40) 
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For  the  Maxwellian  distrlbutlc 


n(E  -  e0^) 


f“En(E) 


=  e  ^ 


vhich  is  identical  to  Eq.  (5),  the 


function. 


result  obtained  from  the  distributic 


The  p„te„tlal  ia  atth.ctlve  re.  p.etena,  ana  an  p.„ro„s  incident  on 
heath  with  initial  energies  between  zero  and  some  critical  energy  E 

satellite.  The  critical  energy  is  that  energy  which  gives  an' 
impact  parameter  of  a  for  a  grazing  collision. 


2  „2 


a  =  R 


-  a 


If  n^(E)  is  the  density  for  protons,  then 
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^  J  ^  ^  dE  (k3) 

R  o  E 

c 

where  the  first  Integral  involves  a  constant  cross  section  times  the  dif¬ 
ferential  flux,  and  the  second  Involves  the  energy-dependent  cross  section 
times  the  same  flux.  The  Integration  may  be  carried  out  for  a  Maxwellian 
distribution  with  the  result  given  by  Eq.  (15).  As  before,  the  result  in 
Eq.  (4-3)  is  true  only  if  the  potential  satisfies  Eq.  (2k), 

If  the  critical  energy  E^  is  much  greater  than  the  characteristic 
plasma  energy  kT,  nearly  all  the  Incident  protons  will  strike  the  satellite. 
Since  from  Eq.  (k2) ,  the  ratio  of  E^  to  kT  is  equal  to  (ix,  the  strong  po¬ 
tential  condition  in  Eq.  (1?)  that  (ix  »  1  follows. 
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VIII.  coNcmsiows 

Most  of  the  difficulties  and  inadequacies  of  ordinary  probe  theory 
are  well  known.  Even  in  the  simplest  theory,  that  of  the  idealized 
sheath,  inconsistencies  are  encountered  in  the  methods  for  the  solution  of 
a  basic  problem;  the  determination  of  the  potential  of  the  floating  probe 
in  a  collisionless,  homogeneous,  isotropic  plasma,  in  this  report  the 
problem  was  treated  from  the  point  of  view  of  distribution  functions,  and 
the  inconsistencies  are  associated  with  the  assumed  boundary  condition 
for  the  Boltzmann  equation  at  the  edge  of  the  sheath.  If  the  inconsistencies 
are  neglected,  the  results  show  that  the  potential  is  almost  independent 
of  the  sheath  thickness. 

Recently,  floating  probe  theories  have  been  adapted  to  the  problems  of 
satellite  charge-up  and  electric  drag.  The  results  of  this  report  show 
that  even  in  the  presence  of  suprathermal  flux  distributions,  such  as  the 

Van  Allen  flux,  the  satellite  potential  is  a  weak  function  of  the  sheath 
thickness. 

Diode  theory  is  a  special  limiting  case  of  idealized  sheath  theory, 
and  its  use  in  sheath  thickness  claculatlons  is  subject  to  considerable 
inaccuracies.  Because  of  the  weak  dependence  of  the  potential  on  the  thick¬ 
ness,  the  use  of  diode  theory  leads  to  little  error  in  satellite  potential 
calculations,  but  its  use  in  more  sensitive  areas,  such  as  in  electric 
drag  claculations,  will  likely  lead  to  significant  errors. 

To  determine  the  satellite  or  probe  equilibrium  potential,  it  is 
convenient  to  calculate  first  the  positive  and  negative  currents  to  the 
probe  as  functions  of  the  probe  potential,  and  then  to  equate  the  currents 
and  solve  for  the  potential.  If  ijr  (e)  is  the 


unperturbed  flux  of  negative 


particles  with  energies  greater  than  E,  then  the  flux  :  of  such  particles 
to  an  absorbing,  negatively  charged  spherical  probe  with  potential  0^^  is 

given  by 


(8) 


where  z  =  -e  ^  ,  and  -e  is  the  charge  on  the  particle.  The  flux  given  by 
Eq.  (8)  is  a  consequence  of  the  idealized  sheath  theory,  and,  although 
it  does  not  depend  on  either  the  probe  or  sheath  dimensions,  it  is  neces¬ 
sary  to  assume  that  the  potential  throughout  the  sheath  region  lies  between 
zero  and  •  However,  if  the  unperturbed  plasma  is  electrically  neutral, 
the  theory  predicts  a  positive  sheath  about  a  negative  probe,  and  the 
condition  is  satisfied.  In  fact,  a  positive  sheath  means  the  potential  is 
a  monotonic  function  of  the  radial  distance  throughout  the  sheath,  and 
lies  between  zero  and  a  limiting  potential  which  varies  as  between 
at  the  inside  edge,  and  zero  at  the  outer  boundary. 

If  the  unperturbed  negative  paritcles  have  a  Maxwellian  distribution, 
Eq.  (8)  reduces  to  the  well-known  expression 

r=  >K0)  exp  j^j 

Another  example,  of  interest  in  connection  with  the  trapped  electrons  in 
the  earth's  radiation  belt,  is  a  flux  i|i(E)  which  goes  as  E  for  E  >  z, 
and  for  which  Eq.  (8)  gives 


_  'l!(z) 

n  +  1 


(9) 


The  positive  particles  are  moving  in  an  attractive  field  and  a  general 
expression  for  the  flux  to  the  probe,  corresponding  to  that  for  the  negative 
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particles  given  by  Eq.  (8),  cannot  be  obtained.  An  upper  limit  to  the 
positive  flux  can  he  obtained  with  the  assumption  that  all  trajectories 
are  filled  for  which  the  particle  energy  and  angular  momentum  at  the  outer 
boundary  of  the  sheath  are  equal  to  the  energy  £ind  angular  momentum  at  the 
inner  edge.  If  the  unperturbed  velocity  distribution  of  positive  particles 
is  not  restricted,  this  assumption  Implies  that  the  potentleO.  tiurou^out 
the  sheath  lies  below  a  limiting  potential  which  varies  as  1  r^  between 
at  the  Inside  edge  and  zero  at  the  outer  boundary.  If  this  condition 
on  the  potential  holds,  and  if  the  positive  particles  have  a  Ifexwellian 
velocity  distribution,  the  positive  flux  to  the  probe  is  given  by  the 
Mott-Smith  expression  shown  in  Eq.  (15). 

If  the  high-energy  electron,  fliix  in  the  more  intense  regions  of  the 

T  2 

earth's  radiation  belt  is  ~  5x10'  electrons/cm  /sec  above  55  kev,  with  the 

_k 

flux  varying  as  E  ,  and  if  the  protons  have  a  Meixwelliein  distribution  at 
~  3000°  K  with  a  density  of  ~  lOO/cm^,  the  satellite  equilibrium  potential 
will,  from  Eqs.  (9)  and  (15),  be  ~  10,000  volts  negative.  This  potential 
should  be  reduced  (in  magnitude)  by  a  factor  of  ~  5  for  a  satellite  with 
a  diameter  much  smaller  than  1  meter,  and  should  be  multiplied  by  ~  3  for 
a  satellite  much  greater  than  1  meter  in  diameter. 
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